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Executive  Summary: 


Described  herein  are  the  key  findings  and  results  associated  with  the  Advance  Microwave 
Ferrite  Research  (AMFeR),  Phase  III  project.  The  purpose  of  this  research  endeavor  is  to 
devise  ferrite  materials  for  microwave,  self-biased  circulator  applications.  Materials 
suitable  for  this  study  are  the  class  of  M-type,  hexagonal  ferrites,  which  typically  have 
large  crystalline  anisotropies  that  tend  align  the  magnetic  moments  of  the  ferrite  in  a 
single  direction  with  the  hope  of  achieving  a  quasi-saturated  state.  In  addition  to  having 
large  crystalline  anisotropy,  the  materials  must  also  possess  large  magnetic  coercivities 
that  keep  the  ferrite  saturated  even  in  the  presence  of  a  strong,  shape-dependent 
demagnetization  field.  As  a  rule  of  thumb,  the  coercivity  should  be  of  the  same  strength 
(or  greater)  as  the  magnetic  saturation  to  achieve  self-bias  operation. 

The  central  task  of  the  project  is  to  fabricate  ferrites  that  have  a  high  magnetic  saturation, 
high  coercivity  and  low  FMR  linewidth.  These  ferrites  are  devised  for  microwave  ferrite 
devices  in  the  context  of  some  integration  process.  Past  reports  for  Phase  I  and  II  explain 
several  techniques  -  including  liquid  phase  epitaxy,  chemical  vapor  deposition  and 
sputtering  -  that  have  been  tried  to  achieve  this  end,  but  have  eventually  been  discarded 
for  various  reasons.  Such  reason  include  slow  deposition  growth  rates,  incompatible 
material  stresses  and  strains,  high  FMR  linewidths,  poor  coercivities,  etc.  The  method 
currently  being  considered  is  based  on  a  ferrite  and  epoxy  compound  that  is  hand  mixed, 
compressed  and  magnetically  aligned  by  external  magnets.  This  method  yields  lower 
magnetic  saturations  (i.e.  2,200  G),  but  sufficient  coercivities  that  keep  the  ferrite  self- 
biased.  Attempts  are  currently  being  made  to  integrate  this  compound  into  a  working 
device.  The  outcome  of  this  activity  will  be  reported  in  the  Phase  IV  Final  Report. 

Concurrent  with  the  material  fabrication  research  is  research  in  the  area  of  wideband, 
high-isolation  circulator  design,  with  particular  attention  being  focused  on  the  integration 
of  a  circulator  with  an  antenna.  The  aim  of  this  research  is  bring  about  an  integrated 
assembly  that  is  light-weight,  low-profile  and  small.  Two  prototype  assemblies  have  been 
designed,  fabricated  and  tested  to  prove  the  soundness  of  the  methodology.  Work  has 
also  been  done  on  the  modeling  of  electromagnetic  waves  in  ferrites  using  a  finite 
difference,  time-domain  methodology.  This  activity  resulted  in  a  high  quality  simulation 
of  a  ferrite  phase  shifter. 

The  key  team  members  of  this  project  are  divided  into  two  functional  groups:  Material 
Science  and  Microwave  Device  Design.  The  Material  Science  Group,  headed  by 
Professors  David  Mcllroy  and  Wei  Jiang  Yeh,  is  responsible  for  ferrite  material 
fabrication,  characterization  and  validation.  The  Microwave  Device  Group,  headed  by 
Professor  Jeffrey  Young,  is  responsible  for  device  design,  fabrication  and  test.  Each 
group  has  prepared  a  report;  the  reports  are  attached  as  appendices.  The  appendices  are 
organized  alphabetically  by  author  as  follows: 

Appendix  A:  Wei  Jiang  Weh  and  David  Mcllroy,  “Fabrication  of  Barium  Ferrite  Disks’7 
Appendix  B:  Jeffrey  L.  Young,  “Microwave  Circulator  Design” 
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In  addition  to  the  key  tasks  associated  with  the  project,  the  Microwave  Device  Group 
continued  its  work  on  the  construction  of  its  2-40  GHz  anechoic  chamber,  1 10  GHz 
probing  station  and  FMR  linewidth  measurement  apparatus.  Each  of  these  projects  was 
completed  in  the  middle  of  Phase  Four  and  will  be  discussed  in  that  report. 

The  work  performed  during  Phase  III  also  led  to  a  partnership  with  Navy  SPA  WAR,  San 
Diego.  Working  with  J.  Rockway,  J.  Allen  and  D.  Arceo,  the  team  devised  an  idea  of 
how  to  use  circulators  in  the  design  of  a  comb  RF  receiver  for  the  Marines.  This  idea  led 
to  the  submission  of  a  joint  patent  application  between  SPAWAR  and  the  University  of 
Idaho.  The  idea  also  led  to  a  three  year  funded  ONR  project  under  the  BAA  07-037. 

Finally,  as  a  matter  of  international  recognition.  Prof.  Jeffrey  Young  was  elected  to  the 
grade  of  IEEE  Fellow  during  Phase  Three  for  his  work  in  wave  propagation  in  complex 
media. 

Future  Work: 

Both  the  material  processes  and  laboratories  are  in  place  to  commence  with  the 
fabrication  of  actual  self-biased  circulators.  During  Phase  Four,  both  groups  will  bring 
their  separate  skills  and  knowledge  together  to  design,  fabricate  and  test  a  microwave, 
self-biased  circulator.  At  the  time  of  this  writing,  the  two  groups  have  obtained  results 
that  look  very  promising.  In  fact,  by  using  strontium,  bulk  ferrite,  the  Microwave  Device 
Group  has  fabricated  a  working  self-bias  circulator.  Further  details  of  these  activities  will 
be  stated  in  the  Phase  Four  Final  Report. 

Journal  Publications: 

R.  S.  Adams,  B.  O’Neil,  and  J.  L.  Young,  “Integration  of  microstrip  circulators  with 
planar  Yagi  antennas  of  several  directors,”  IEEE  Transactions  on  Antennas  and 
Propagation  (accepted  for  publication;  see  attachments). 

R.  S.  Adams,  B.  O’Neil  and  J.  L.  Young,  “The  circulator  and  antenna  as  a  single 
integrated  system,”  IEEE  Antennas  and  Wireless  Propagation  Letters  (accepted  for 
publication;  see  attachments). 

R.  S.  Adams,  “A  generalized  time-stepping  algorithm  for  electromagnetic  propagation  in 
ferrite  media,”  IEEE  Transactions  on  Antennas  and  Propagation  (accepted  for 
publication;  see  attachments). 
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Conference  presentations: 


W.  J.  Yeh,  A.  Abuzir  and  Y.  Kranov,  “Barium  hexaferrite  films  made  by  sputtering  and 
liquid  phase  epitaxy  method  for  microwave  application,”  Sixteenth  International 
Conference  on  Composite  and  NanoEngineering,  July  2007,  Haikou,  China. 

W.  J.  Yeh,  C.  Blengeri-Oyarce,  S.  Pillamari,  J.  M.  Appikonda,  L.  Diaz,  Y.  Kranov  and 
D.  Mcllroy,  “Fabrication  of  barium  ferrite  thick  films,”  l(fh  Annual  Meeting  of  Northwest 
Section  of  American  Physical  Society,  May  2008,  Portland,  Oregon. 

S.  Chava,  W.  J.  Yeh  and  J.  Tilles,  “Fabrication  of  metallic  nanoporous  films,”  American 
Physica  Society  March  Meeting,  March  2008,  New  Orleans,  Louisiana. 

W.  J.  Yeh,  S.  Chava,  J.  Tilles,  “Metallic  nanoporous  films  fabricated  by  etching,”  9th 
Annual  meeting  of  Northwest  Section  of  American  Physical  Society,  May  2007,  Pocatello, 
Idaho. 

W  J.  Yeh,  S.  Chava,  J.  Tilles,  “Gold  and  lead  metallic  nanoporous  films  made  by 
etching,”  at  American  Physical  Society  March  Meeting,  March  2007,  Denver,  Colorado. 

R.S.  Adams,  J.L.  Young  and  B.  O’Neil,  “Design  of  microstrip  antenna/ferrite  circulator 
systems,”  IEEE  Antennas  and  Propagation  International  Symposium,  Honolulu,  Hawaii, 
June  2007. 

R.S.  Adams  and  J  L  Young,  “The  finite  difference  time  domain  algorithm  applied  to 
ferrite  media,”  IEEE  Antennas  and  Propagation  International  Symposium,  Honolulu, 
Hawaii,  June  2007. 

Theses  and  Dissertations: 

Benton  O’Neil,  “Self-biased  microstrip  Y-junction  circulator  design  using  effective 
cavity  models  and  power  ratios,”  Department  of  Electrical  and  Computer  Engineering, 
University  of  Idaho,  Moscow,  Idaho,  July  2007. 

Patents: 

J.  W.  Rockway,  J.  C.  Allen,  D.  Arceo  and  J.  L.  Young,  “Broadband  Antenna  Comb 
Circulator  CLC/CLAC  (BA-C4)”  (Submitted  July  2007). 
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Appendix  A 

Fabrication  of  Barium  Ferrite  Disks 

Wei  Jiang  Yeh  and  David  N.  Mcllroy 


1.  Abstract 

During  recent  years  the  need  for  high  quality  self-biased  barium  ferrite  (BaFe120„)  thick 
films  had  been  increasing  due  to  their  chemical  stability  and  oriented  hexagonal 
crystalline  structure.  Our  goal  is  to  fabricate  barium  ferrite  thick  disks  to  be  incorporated 
in  self-biased  microwave  devices.  Different  methods  such  as  sputtering,  pulse  laser 
deposition,  chemical  vapor  deposition  and  modified  liquid  phase  deposition  have  been 
used  with  limited  to  no  success  for  thickness  above  300pm.  Excessive  residual  stresses  of 
deposited  BaM  films  films  and/or  low  coercivities  are  the  main  problems  to  be  solved  by 
this  project.  We  present  a  low  cost  solution  using  BaFe,jO„  nanopowder  mixed  with 
epoxy.  The  mix  is  placed  on  alumina  substrates  to  fabricate  500pm  thick  disks  of  BaM  to 
achieve  good  magnetic  properties  such  as  magnetic  saturations  from  2,000  to  2,500  G 
and  a  coercivities  fom  3800  to  4000  Oe.  In  addition,  we  have  integrated  these  barium 
ferrite  thick  disks  into  microwave  devices  and  successfully  deposited  copper  contact  lines 
for  testing  and  use. 

II.  Proposal  Objective 

•  Develop  barium  ferrite  thick  disks  to  be  incorporated  in  self-biased  microwave 
devices. 

•  Integrate  the  barium  ferrite  thick  disk  within  the  substrate  and  pattern  the  metallic 
traces  to  achieve  a  microwave  circuit. 

III.  Methodology 

Fabrication  of  the  Barium  Ferrite  thick  film 

Each  barium  ferrite  thick  disk  is  made  from  a  combination  of  commercially  available 
barium  ferrite  nanopowder  (barium  dodecairon  nonadecaoxide  BaFeiA?)  ar,d  30  minutes 
epoxy  with  a  ratio  of  80/20,  respectively.  The  barium  ferrite  nanopowder  is  weighed  and 
placed  into  the  mixing  container;  the  epoxy  is  then  weighed  and  added  to  the  mix.  To 
help  with  the  mixing  process  1  ml  of  acetone  is  added  to  the  mix.  This  first  part  of  the 
mixing  is  done  for  2  minutes  and  the  goal  is  to  incorporate  all  the  material  together  to 
form  a  dough-like  compound,  as  shown  in  Figure  1.  After  that,  the  dough  is  mixed  by 
hand  for  another  2  minutes  to  achieve  a  uniform,  homogeneous  composition;  see  Figure 

2. 
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A  small  piece  of  dough  is  placed  in  a  mechanical  press  containing  one  inch  magnets.  The 
mixed  material  is  placed  in  the  substrate  and  carefully  pressed,  as  shown  in  Figure  3;  then 
another  piece  of  the  mixed  substance  is  added  on  top.  The  press  is  closed  carefully,  as 
shown  in  Figure  4.  Each  press  contains  two,  one  inch  magnets,  with  one  on  top  and  one 
on  bottom.  The  function  of  the  magnets  is  to  align  the  magnetic  particles  of  the  mixed 
substance  so  that  their  c-axes  are  uniformly  parallel  and  oriented.  The  samples  are  kept 
in  the  high  magnetic  field  for  48  hours,  and  then  released  and  cured  at  room  temperature 
for  another  48  hours. 


Polishing 

After  the  curing  process  is  completed,  each  sample  is  polished  using  a  precision 
lapping/polishing  machine  (UNIPOL-810),  as  shown  in  Figure  5.  The  polishing  process 
is  done  in  several  steps  with  the  lapping  machine  spinning  at  250  rpm: 

1.  Samples  are  polished  with  320  waterproof  sandpaper  to  remove  the  excess 
material. 

2.  After  all  the  excess  material  is  removed,  a  polishing  cloth  is  placed  in  the 
polishing  machine.  Three  different  polishing  slurries  are  sequentially  used  -  3.0 
micron,  0.5  micron  and  a  glass  polish  0.05  micron.  Each  of  the  slurries  is  used  for 
6  minutes  for  each  side. 
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3.  After  the  polishing  is  done  the  samples  are  place  on  a  hot  plate  at  45°C  for  24  hrs 
to  dry. 


Figure  5.  Sample  being  polish 


Gap  Filling 

During  the  drying  step  some  gaps  develop  between  the  ferrite  disk  and  substrate.  To 
remedy  this  problem  an  epoxy  filler  is  prepared  and  placed  into  the  gaps  using  a  thin 
blade;  this  particular  step  requires  extreme  care.  It  is  important  that  the  epoxy  be  not 
applied  inside  the  via  holes  that  are  used  as  part  of  the  electrical  circuit.  Once  all  the 
excess  material  is  removed,  the  samples  are  placed  again  on  the  hot  plate  for  24  hrs  for 
curing. 

Cleaning  and  Spin  Coating  Process 

Once  the  samples  are  cured,  each  one  is  then  cleaned  with  acetone  and  alcohol  using 
Kimwipes;  they  are  subsequently  dried  with  compressed  air.  To  achieve  a  uniform 
surface  for  copper  patterning,  each  sample  is  placed  in  the  spin  coater  K.W-4A,  shown  in 
Figure  6,  and  coated  with  a  UV  adhesive  (Norland  Optical  Adhesive  68).  The  spin 
process  lasts  about  one  minute  at  a  rate  4,000  rpm.  Samples  are  then  exposed  in  a 
30mW/cm2  UV  light  for  15  minutes  and  placed  on  a  hot  plate  at  50°C  for  12  hrs.  The 
same  process  is  performed  on  both  sides  of  the  substrate. 
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Sputtering  Process 


Figure  6  Spin  Coater  KW-4A 


Sputtering  is  a  technology  in  which  a  material  is  released  from  a  pure  material  source 
named  target,  such  Au,  Cu,  or  Ti.  The  substrate  is  placed  in  a  vacuum  chamber  with  the 
target  and  an  inert  gas,  such  as  argon,  is  introduced  at  low  pressure,  typically  in  the  range 
of  10-100mTorr.  Gas  plasma  is  created  by  using  a  power  source,  which  causes  the  gas  to 
ionize.  The  ions  are  accelerated  towards  the  surface  of  the  target  causing  atoms  of  the 
source  material  to  break  off  from  the  target,  which  then  vaporize  and  condense  on  all 
surfaces  including  the  substrate.  Difference  in  pressure  and  voltage  typically  relate  to  the 
manor  in  which  the  ion  bombardment  of  the  target  is  realized.  A  schematic  diagram  of  a 
typical  DC  sputtering  system  is  shown  in  Figure  7. 


Sputtering  Target  (gun) 


Substrate  Holder 

Figue  7.  Graphic  representation  of  the  sputtering  process 


We  used  a  single  gun  sputtering  system.  Argon  DC  plasma  deposition  is  used  to  fabricate 
all  the  layers  of  titanium  and  copper.  The  deposition  was  done  at  a  pressure  of  8mTorr 
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with  a  DC  power  input  of  150W.  At  these  settings,  the  deposition  rate  was  0.26pm/min 
for  copper  and  500A/min  for  titanium. 


Figure  8.  Sample  after  being  sputtered 


Measuring  Samples 

The  ferrite  samples  are  characterized  using  a  vibrating  sample  magnetometer  (VSM),  as 
shown  in  Figure  9.  Two  measurements  are  taken:  one  at  0-degree  angle  (out  of  plane)  and 
at  90-degree  angle  (in  plane).  Because  the  hexagonal  substrates  cannot  be  placed  on  the 
VSM  to  be  measured,  circular  samples  are  made  at  the  same  time  from  the  same  ferrite 
compound  as  the  original  samples.  This  is  shown  in  Figure  10. 


Figure  9.  Vibrating  Sample  Magnetometer 
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Circular  sample  unpolished  Circular  sample  polished 


Figure  10.  Circular  samples  to  be  measure  by  the  VSM 

IV.  Results 

Over  the  last  year  different  nanopowder/epoxy  ratios  have  been  tested.  As  it  can  be 
observed  in  Graphic  1,  the  80/20  ratio  has  the  best  overall  results.  A  90/10  ratio  was  also 
tested,  but  the  mixing  process  was  difficult  and  the  uniformity  of  the  compound  was 
compromised.  With  the  80/20-recipe  ratio  we  obtained  very  manageable  and  uniformed 
material,  which  resulted  in  similar  and  consistent  magnetic  properties  for  all  samples. 


Overall  Ratio  Effectiveness  on  VSM  Parameters 


Powder  Ratio 


■  COE  ■4piMS  □  MR  BSQ 

Graphic  1 .  Recipe  ratios  tested  along  with  overall  outcome  of  all  measure  parameters 


Different  puck  thicknesses  have  also  been  tested  from  0.40  mm  to  almost  3.00  mm  thick. 
Graphic  2  shows  that  thicker  barium  ferrite  films  tend  to  have  better  overall  performance. 
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Our  objective  was  to  have  barium  ferrite  disks  of  0.5  mm  in  thickness  with  good 
magnetic  properties.  The  data  confirms  that  this  has  been  accomplished. 


Magnetic  Porperties  vs.  Thickness 


0.49 


0.82 


1.07  1.53  2.02 

Thickness  (mm) 


2.56 


2.92 


I  COE  B4piMS  DMR  DSQ 


Graphic  2.  Different  puck  thickness  and  its  magnetic  properties  performance. 


Magnetic  Properties  Measure  bv  the  Vibrating  Sample  Magnetometer  (VSM) 

Typical  measured  data  obtained  by  the  vibrating  sample  magnetometer  is  shown  in 
Figure  11.  The  most  important  magnetic  properties  that  we  take  into  consideration  for 
this  project  are  the  squareness  of  the  hysteresis  loop;  the  second  property  is  the  coercive 
field  that  indicates  the  strength  of  the  magnetic  field  needed  to  reduce  the  magnetization 
to  zero.  The  third  property  is  the  saturation  magnetization  or  4jiMs.  which  indicates  the 
maximum  possible  magnetization  of  the  sample.  Lastly  the  forth  property  is  the  remanent 
magnetization,  which  indicates  the  magnetism  that  remains  after  a  saturating  applied  field 
is  removed. 
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In  the  data  shown  below,  we  can  observed  that  the  samples  presents  a  squareness  of  0.89 
S.  a  coercive  field  of  4127  Oe.  a  4jiMs  of  2.469  Gauss  and  a  remanent  magnetization  of 
2.201  Gauss.  For  self-biased  operation,  the  coercivity  should  exceed  the  saturation,  which 
we  see  is  the  case.  Thus,  this  sample  is  a  good  candidate  for  a  self-biased  circulator. 


Figure.  1 1 .  VSM  measurements  for  the  in-plane  (red  line)  and  out-of-plane  (blue  line) 
hysteresis  loops  for  barium  ferrite  thick  film. 


Table  1  shows  the  results  of  200  samples  measured  by  the  VSM.  Judging  by  the  standard 
deviation,  we  can  conclude  that  the  technique  described  herein  has  high  rate  of 
repeatability  and  consistently  produces  good  thick  disks  with  a  squareness  average  of 
0.83  and  a  coercivity  average  of  3,986  Oe.  Samples  also  show  good  values  for  the 
saturation  magnetization  with  an  average  of  2,255  G  and  1,862  G  for  the  remanent 
magnetization. 


Thickness 

(mm) 

Squamess 

(Sq) 

Coercivity 

Field 

Saturation 
Magnetization 
(4  n  Ms) 

Remanent 

Magnitization 

(MR) 

Mean  (p) 

0.76 

0.83 

3986 

2,255 

1,862 

STDV  («) 

0.24 

0.03 

296 

107 

81 

Table  1 .  Value  of  the  mean  and  standard  deviation  of  200  samples  measured  by  the  VSM 
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Appendix  B 


Microwave  Circulator  Design 

Jeffrey  L.  Young 


I.  Abstract 

During  Phase  III,  a  significant  portion  of  the  Group’s  energy  was  spent  on  developing 
techniques  for  the  integration  of  a  microstrip  circulator  with  a  planar  Yagi  antenna  to 
create  a  compact,  light-weight,  low-profile  integrated  subassembly.  As  reported  here,  this 
research  activity  was  very  successful  with  the  outcome  of  building  two  such,  fully- 
functioning  assemblies.  Additionally,  the  team  also  developed  a  numerical  algorithm  that 
predicts  the  time-domain  wave  behavior  of  a  microwave  signal  within  a  ferrite.  The 
device  that  was  used  for  verification  of  the  algorithm  was  a  ferrite  phase  shifter  that 
exhibited  excellent  return  loss  and  linear  phase  characteristics. 

II.  Proposal  Objectives: 

The  primary  objectives  for  Phase  III  were: 

•  To  develop  a  wideband,  high  isolation,  self-biased  circulator  for  microwave  signal 
applications; 

•  To  develop  tools  and  models  for  the  design  and  development  of  ferrite  devices; 

•  To  integrate  ferrite  devices  with  other  microwave  devices  (e.g.  antennas)  to  form 
a  light-weight,  low-profile  subassembly. 

As  discussed  herein,  the  last  two  objectives  have  been  fully  realized.  As  for  the  first 
objective,  we  did  not  achieve  full  success  until  the  middle  of  the  fourth  phase  of  the 
project;  this  will  be  reported  in  the  Phase  Four  Final  Report.  One  of  many  reasons  for  this 
delay  was  vender-related.  A  vender  that  was  contracted  to  fabricate  a  specialized 
strontium  ferrite  and  other  venders  who  did  tasks  like  polishing  and  etching  turn  what 
was  thought  to  be  a  six  week  process  into  a  six  month  process. 

III.  Methodology: 

The  methodologies  and  results  of  the  Microwave  Device  Group  are  described  in  the 
following  papers: 

•  R.  S.  Adams,  B.  O’Neil,  and  J.  L.  Young,  “Integration  of  microstrip  circulators 
with  planar  Yagi  antennas  of  several  directors.” 

•  R.  S.  Adams,  B.  O’Neil  and  J.  L.  Young,  “The  circulator  and  antenna  as  a  single 
integrated  system.” 

•  R.  S.  Adams,  “A  generalized  time-stepping  algorithm  for  electromagnetic 
propagation  in  ferrite  media.” 
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Note:  These  papers  area  appended  to  this  report.  All  three  of  these  papers  have  been 
accepted  for  publication.  Unfortunately,  attribution  to  the  Phase  One  contract  number 
was  incorrectly  cited  by  the  lead  author;  the  research  associated  with  these  papers 
commenced  towards  the  end  of  Phase  Two  (but  not  reported  in  the  Phase  Two  Final 
Report)  and  was  carried  out  throughout  Phase  Three. 

IV.  Results 

See  the  aforementioned  citations  and  comments. 

V.  Future  Work 

Phase  Four  is  the  final  phase  of  the  Advanced  Microwave  Ferrite  Research  (AMFeR) 
project.  As  mentioned  in  the  Executive  Summary,  a  self-biased  circulator  has  been 
fabricated  and  tested.  The  results  of  those  tests  clearly  indicate  that  reasonable  insertion 
losses  on  the  order  of  1.5  dB  and  isolations  on  the  order  of  15  dB  are  possible  over  a 
narrow  band  when  the  device  is  operating  at  22  GHz.  The  goal  of  Phase  Four  is  to 
improve  the  bandwidths  without  compromising  the  insertion  and  return  loss  values. 
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Integration  of  a  Microstrip  Circulator  with  Planar 
Yagi  Antennas  of  Several  Directors 

Ryan  S.  Adams,  Member,  IEEE,  Benton  O’Neil,  Student  Member,  IEEE  Jeffrey  L.  Young,  Senior  Member,  IEEE 


Abstract—  Two  Yagi  planar-type  antennas  are  presented  that 
exhibit  increased  gain  and  decreased  beamwidth  in  the  principle 
planes  when  compared  with  the  baseline  antenna  from  the  open 
literature.  Two  systems  are  also  presented  in  which  each  planar 
antenna  is  integrated  with  a  ferrite  circulator.  Excellent  pass  band 
characteristics  are  observed  in  the  experimental  and  numerical 
data  for  both  systems.  Further,  one  of  the  systems  is  magnetically 
tuned  to  provide  a  superb  30  dB  isolation  over  a  bandwidth 
of  850  MHz.  The  other  system  is  also  magnetically  tuned  to 
provide  30  dB  isolation  over  a  300  MHz  bandwidth  and  40  dB 
isolation  over  a  100  MHz  bandwidth.  Both  systems  maintain  an 
excellent  return  loss  specification  and  show  improved  insertion 
loss  specifications  even  though  the  ferrite  is  tuned  to  be  in  partial 
saturation.  Extensive  measurement  data  of  the  ferrite  material 
and  system  network  parameters  are  provided  to  corroborate  the 
claims  made  herein. 

Index  Terms —  Ferrites,  Microwave  Circulators,  Microstrip 
Antennas. 


Introduction 

A  typical  monostatic  RADAR  system  consists  of  a  trans¬ 
mitter,  receiver,  antenna  and  a  circulator,  as  shown  in  Figure 
1.  The  design  of  these  systems  usually  consists  of  separating 
the  circulator  from  the  remaining  components  and  designing 
it  to  some  standard  interface  specification  (e.g.,  50D),  as  is 
the  adjoining  circuitry.  Each  component  is  then  combined  to 
form  the  system.  Although  this  approach  will  yield  adequate 
results,  other  approaches  can  provide  greater  performance  and 
flexibility.  We  present  herein  a  methodology  that  considers  the 
design  of  a  planar  antenna  with  a  microstrip  circulator  as  a  sin¬ 
gle,  integrated  component.  The  antenna  topology  considered 
is  the  planar  Yagi-Uda  array  [1].  To  validate  the  proposed 
methodology,  two  integrated  antenna/circulator  systems  are 
presented. 

Simulation  data  are  provided  for  three  substrate  antennas, 
baseline,  two  director  and  three  director  antennas.  The  baseline 
antenna  is  derived  from  the  Yagi-Uda  antenna  presented  by 
Kaneda  et.  al.  [1].  This  antenna  exhibits  a  fairly  narrow 
beamwidth  in  the  E- plane,  but  the  //-plane  beamwidth  is 
very  wide.  The  performance  of  this  baseline  antenna  topology 
is  enhanced  through  the  use  of  additional  director  elements  that 
provide  a  significant  improvement  in  //-plane  beamwidth.  The 
directivity  of  the  antenna  is  also  improved  as  a  consequence 
of  the  modifications  presented  herein.  When  these  systems 
are  utilized  for  communications  or  RADAR  applications,  the 

Manuscript  received  XXXX  XX,  20XX;  revised  XXXX  XX,  20XX.  This 
work  was  supported  by  the  Office  of  Naval  Research  under  Award  Number 
N000 14—04— 1-0272. 

The  authors  are  with  the  Department  of  Electrical  and  Computer  Engineer¬ 
ing,  University  of  Idaho,  Moscow  Idaho  (email:  radams@uidaho.edu). 


Fig.  1.  A  monostatic  RADAR  system. 


improved  gain  compensates  for  the  inherent  insertion  loss 
associated  with  the  circulator. 

Baseline  Antenna 

We  use  as  a  baseline  a  modified  form  of  the  antenna  of 
Kaneda;  modifications  were  made  to  account  for  our  choice 
of  substrate  thickness  (0.5mm)  and  substrate  permittivity  (4.5). 
This  antenna  topology  is  shown  in  Figure  2  and  has  the  fol¬ 
lowing  geometrical  parameters:  Ssub  =  1.2mm,  La r  —  3.3mm, 
War  =  0.6mm  and  Sd,r  =  2.4mm.  All  other  geometrical 
parameters  are  identical  to  that  of  Kaneda,  except  for  minor 
changes  in  the  miters.  This  antenna  was  simulated  to  determine 
pattern  and  impedance  information;  Figures  3  and  4  show 
the  E-plane  and  //-plane  patterns  respectively  and  Figure  5 
shows  the  input  impedance  to  the  antenna.  The  patterns  exhibit 
the  following  performance:  E-plane  beamwidth  is  87°  with 
sidelobe  levels  of  16.5  dB;  antenna  gain  is  5.12  dB;  //-plane 
beamwidth  is  greater  than  180°  with  no  sidelobes.  The  real 
part  of  the  input  impedance  undulates  between  25  and  50  fl  in 
the  12  to  17  GHz  band.  The  imaginary  part  undulates  between 
20  and  -20  fi  in  the  same  band.  Provided  that  the  circulator 
presents  to  the  antenna  an  impedance  that  is  close  to  these 
values,  a  good  impedance  match  between  the  circulator  and 
the  antenna  can  be  designed. 

This  antenna  exhibits  adequate  performance  in  the  Ku- 
band;  however,  the  //-plane  beamwidth  is  quite  broad.  To 
allow  tradeoffs  in  system  design,  we  modify  the  topology  to 
include  additional  director  elements,  which  improves  //-plane 
beamwidth  as  well  as  antenna  gain,  as  described  next. 


Double  Director  Antenna 

Consider  the  effect  of  an  additional  director  element  to 
the  geometry  of  the  baseline  antenna,  as  shown  in  Figure 
6.  To  be  consistent  with  the  baseline  antenna,  we  chose  the 
following  modified  geometrical  parameters:  t  =  1.2mm, 

Ldir  =  3.3mm,  War  -  0.6mm,  Sdir  =  2.4mm  and  Sdir-dir 
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Fig.  2.  Single  director  planar  Yagi-LJda  antenna. 


Fig.  5.  Input  impedance  of  the  baseline  antenna. 


Fig.  6.  Double  director  planar  Yagi-Uda  antenna. 


Fig.  3.  E-plane  pattern  comparison  for  the  baseline,  double  director  and 
triple  director  antennas. 


Fig.  4.  H -plane  pattern  comparison  for  the  baseline,  double  director  and 
triple  director  antennas. 


=  2.4mm.  The  remaining  parameters  are  identical  to  those 
of  the  baseline  antenna.  The  E-plane  and  //-plane  antenna 
patterns  were  simulated  with  results  shown  in  Figures  3  and  4. 
From  these  plots,  we  note  the  following:  E-plane  beamwidth 
is  93°  with  sidelobe  levels  of  13.9  dB;  antenna  gain  is  5.86 
dB;  //-plane  beamwidth  is  140°  with  no  sidelobes.  The  input 
impedance  to  this  antenna  was  approximately  the  same  as  that 
of  the  baseline  antenna,  and  hence  the  simulated  impedance 
is  not  shown. 

Triple  Director  Antenna 
The  final  antenna  under  consideration  uses  a  third  director 
element  to  the  geometry  of  Figure  6.  Again,  we  maintain 
consistency  with  the  previous  antennas  by  choosing  the  fol¬ 
lowing  modified  geometrical  parameters:  SSub  =  1.2mm,  Ldir 
=  3.3mm,  Wdir  =  0.6mm,  Sdir  —  2.4mm  and  S  dir -dir  ~ 
2.4mm.  The  E-plane  and  //-plane  patterns  of  this  antenna 
were  also  simulated  and  the  data  are  shown  in  Figures  3 
and  4.  The  following  features  are  noted  from  these  plots: 
E-plane  beamwidth  is  76°  with  sidelobe  levels  of  13.58  dB; 
antenna  gain  is  7.05  dB;  //-plane  beamwidth  is  105°  with  no 
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sidelobes.  The  input  impedance  of  this  antenna  was  also  very 
similar  to  that  of  the  baseline  antenna,  and  is  not  shown.  This 
topology  yields  best  overall  gain,  and  narrowest  beamwidth  in 
both  principle  planes,  at  the  expense  of  slightly  higher  sidelobe 
levels. 


Antenna  Comparison 

With  the  pattern  data  of  Figures  3  and  4  now  known,  we 
have  several  tradeoffs  to  consider.  First,  the  baseline  antenna  is 
the  most  compact  and  has  the  lowest  sidelobe  levels;  however, 
it  has  lower  gain  and  broader  //-plane  beamwidth.  The  double 
and  triple  director  antennas  each  exhibit  improved  antenna 
gain  and  much  narrower  beamwidth  in  the  //-plane  at  the 
expense  of  higher  sidelobe  levels  in  the  /£-plane  as  well  as 
a  larger  footprint.  These  tradeoffs  can  be  exploited  in  overall 
system  design.  For  example,  the  insertion  loss  (which  is  due 
to  material  loss  in  the  ferrite)  of  a  circulator  is  typically  on  the 
order  of  0.5-1  dB.  An  improvement  in  antenna  gain  by  greater 
than  1  dB  can  offset  the  negative  effects  of  the  circulator  in 
the  overall  system  design. 


Circulator  Design 


To  integrate  a  circulator  with  the  antennas  presented  pre¬ 
viously,  the  ferrite  circuit  needs  to  be  described  in  terms 
of  some  network  parameters.  We  do  this  by  appealing  to 
the  approach  of  Bosma  [2],  who  provided  the  approximate 
impedance  parameters  of  the  ferrite  circuit.  With  this  approach, 
the  ferrite  region  is  assumed  to  be  a  homogeneous  closed 
cavity  with  mixed  PEC  and  PMC  boundary  conditions,  and 
hence  a  trans— impedance  Green’s  function  can  be  defined  as 
follows  [3]: 


G(4>,<t>') 


_  Z/Ioilfa) 

~  27i7£(7/a) 

Z,  /n(7/«)^n(^^') 

2tt  ^  /;(7/a)  +  ("*/M/)f»(7/«)/(7/«) 


Z{  ^ 

+  2tt  ^  -  (wt/M/)^n(7/a)/(7/«)’ 


(1) 


The  Green’s  function  relates  the  vertical  electric  field  on  the 
boundary  of  the  the  ferrite  to  the  azimuthal  magnetic  field 
on  the  walls  of  the  ferrite  cavity.  The  appropriate  impedance 
parameters  can  be  determined  from  these  fields  with  equations 
similar  to 

Z\\  =  /  G(-(  tt/3),0W.  (2) 

/(-(ir/3)-^) 

Although  this  approach  is  well-known  and  extensively  used 
to  design  circulators,  certain  errors  still  exist  in  the  solution. 
To  obtain  a  more  accurate  set  of  network  parameters,  the 
geometry  must  all  be  simulated  via  a  full-wave  electromag¬ 
netic  solver.  However,  full-wave  solvers  take  a  very  long 
time  to  arrive  at  the  desired  solution  when  compared  with 
the  time  required  for  the  Green’s  function.  Hence,  a  two  step 
process  is  required:  first,  the  Green’s  function  is  repeatedly 
invoked  to  find  a  solution  that  is  approximately  correct  (say, 
20%),  and  second  a  full-wave  solver  is  used  to  validate  that 


Fig.  7.  Traditional  microstrip  impedance  equalizer  using  unbalanced  stubs. 


Fig.  8.  Balanced  microstrip  impedance  equalizer 


the  ferrite  circuit  is  adequate,  and  to  generate  more  accurate 
network  parameters  for  use  in  the  system  integration.  Once  the 
network  parameters  are  accurately  known,  the  ferrite  region 
can  be  regarded  as  a  lossless,  nonreciprocal,  three-port,  linear 
network  (LNTPLN)  as  presented  in  [4].  With  this  network 
theory  representation,  the  power  flow  through  the  network  can 
be  evaluated  and  metrics  such  as  return  loss,  isolation  and 
insertion  loss  at  each  port  can  be  determined.  This  is  described 
next. 


System  Integration 

To  create  the  integrated  components,  complcx-to-complex 
impedance  equalizers  need  to  be  designed  to  provide  the  inter¬ 
face  between  each  component  of  the  system  to  optimize  return 
loss,  isolation  and  insertion  loss.  This  approach  of  designing 
the  overall  system  potentially  allows  for  better  performance 
and  flexibility  than  the  common  approach  of  isolating  each 
component  from  the  rest  of  the  system  and  designing  each 
port  to  match  to  some  standard  characteristic  impedance,  such 
as  50Q. 

The  traditional  approach  to  designing  microstrip  impedance 
equalizers  consists  of  a  standard  topology  of  alternating 
through  traces  and  individual  stubs,  as  shown  in  Figure  7.  For 
most  applications,  this  approach  is  sufficient.  However,  when 
driving  an  antenna  like  the  ones  described  previously  with  an 
equalizer  like  that  of  Figure  7,  we  note  that  a  field  imbalance 
will  occur  in  the  transverse  plane  of  the  microstrip  as  shown 
in  Figure  9.  When  this  signal  is  presented  to  the  antenna,  the 
balun  will  be  less  effective  and  hence  the  radiation  and  input 
impedance  will  be  affected.  Consequently,  we  have  chosen 
to  design  our  impedance  equalizers  with  electrically  balanced 
stubs  as  shown  in  Figure  8. 
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Fig.  9.  The  field  structure  m  the  vicinity  of  an  unbalanced  stub. 


Fig.  10.  A  block  diagram  of  the  antenna/circulator  system. 


To  design  the  equalizer  networks,  we  follow  the  approach 
outlined  in  [4]  with  a  few  modifications.  Rather  than  use  the 
transducer  power  gain  as  a  metric  in  conjunction  with  the 
circulator  impedance,  we  choose  to  compute  the  power  flow 
through  the  network.  For  example,  we  know  that  two  of  the 
equalizers  are  terminated  with  a  50f2,  purely  real  impedance, 
and  the  third  equalizer  is  terminated  with  the  impedance 
function  shown  in  Figure  5.  We  choose  a  set  of  lengths  and 
widths  for  each  of  the  elements  shown  in  Figure  8  and  compute 
Ti,  T2  and  T3,  which  are  the  reflection  coefficients  that  load 
the  circulator  at  ports  1,  2  and  3,  respectively;  see  Figure 
10.  Using  these  reflection  coefficients,  we  can  compute  the 
insertion  loss  at  port  1  according  to 


Rh  = 


(i-|r1|2)(i-|rin,1|2) 


|i-ririB>1|2 

the  insertion  loss  between  ports  3  and  1  according  to 

2  Zoi  (i  -  |r3|2)(i  -  in  |2) 

Zo3  |i-rxrtatlp 


il3  i  = 


K-T 


and  the  isolation  between  ports  2  and  1  according  to 

2  r*  lr^  |2\/1  _  ip.  |2\ 


IS2 1  = 


Zoi  (l-IWQ-ir.i2) 

Z02  ii-nr^i2  ’ 


(3) 


(4) 


(5) 


where  I\n>i  is  the  input  reflection  coefficient  at  port  1  of  the 
circulator  when  ports  two  and  three  are  loaded  with  T 2  and  T 3 
respectively.  Here  Zci,  Zo2  and  Z0 3  are  the  port  characteristic 
impedances.  We  can  compute  the  return  loss,  insertion  loss 
and  isolation  for  each  of  the  remaining  two  ports  in  a  similar 


manner.  We  then  compare  the  computed  return  loss,  insertion 
loss  and  isolation  with  the  specificied  acceptable  values.  If 
the  computed  values  satisfy  the  specification  over  a  broader 
bandwidth  than  any  previous  topology,  we  keep  the  design;  if 
not  we  choose  a  new  set  of  lengths  and  widths  for  the  topology 
of  Figure  8  and  reiterate. 


Results 

Two  antenna/circulator  systems  were  designed  according  to 
the  procedure  outlined  above.  Each  design  utilized  a  different 
antenna  element  as  shown  in  Figures  11  and  12.  For  each 
system,  impedance  equalizers  were  required  for  the  two  ports 
not  interfacing  with  the  antenna  to  convert  from  a  standard  50 
A  termination  to  an  impedance  that  optimizes  metrics  of  return 
loss,  isolation  and  insertion  loss.  For  the  circulator/antenna 
port,  the  impedance  equalizer  is  one  of  a  complex-to-complex 
impedance  match,  where  the  input  impedance  of  each  antenna 
is  modified  to  optimize  these  same  metrics.  We  stress  that 
this  approach  is  fundamentally  different  from  designing  the 
circulator  and  antenna  as  separate  50  Q  devices. 


Fig.  1 1.  C ircu lator/an tcnna  system  with  two  director  elements. 


Fig.  12.  Circulator/antenna  system  with  three  director  elements 
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The  two  director  system  was  fabricated  on  0.5mm  thick 
Trans-Tech  D4  material,  which  has  a  dielectric  constant  of  4.5. 
The  ferrite  material  selected  was  Trans-Tech  TTI-3000,  which 
has  the  following  properties:  47 xMs  =  3000  G,  6/  =  12.9, 
tan  6  —  .0005  and  radius  =  1.75mm.  The  externally  applied 
field  required  for  these  material  properties  and  geometry  to 
assure  a  zero  internal  field  was  set  at  2575  Oe.  External  biasing 
for  this  material  was  provided  by  a  laboratory  electromagnet. 

The  externally  measured  hysteresis  loop  for  this  material  is 
shown  in  Figure  13,  where  the  “as-designed”  bias  point  of 
2575  Oe  is  shown  at  the  edge  of  saturation.  The  peak-to- 


s 

i 


Externally  Applied  Field  (Oe) 


Fig.  1 3.  Measured  hysteresis  curve  of  TT 1  -3000. 

peak  linewidth  of  the  ferrite  material  was  also  measured  with 
results  shown  in  Figure  14.  The  peak  to  peak  linewidth,  at  10 
GHz,  is  140  Oe,  which  corresponds  to  a  3  dB  linewidth  of 
238  Oe.  The  published  linewidth  from  Trans-Tech  is  228  Oe 
as  measured  at  9.4  GHz. 


Fig.  14.  Measured  pcak-to-peak  linewidth  of  TT  1-3000  at  10  GHz. 

The  impedance  equalizers  connected  to  the  two  50Q  ports 
of  the  two  director  system  are  identical  and  have  the  following 
geometrical  parameters:  W\  =  0.94mm,  W2  =  0.216mm,  L2  = 
1.142mm,  W3  =  1.097mm,  L3  =  1.067mm,  W4  =  0.216mm, 
L4  =  4.846mm,  Wb  =  2.255mm  and  L5  =  4.219mm.  The 
matching  network  connecting  the  circulator  to  the  two-director 
antenna  has  the  following  geometrical  parameters:  W\  - 


1.2mm,  Li  =  1.201mm,  W2  =  0.216mm,  L2  =  1.214mm,  Wz 
=  2.196mm,  L3  =  2.372mm,  W4  =  0.216mm,  L4  =  6.983mm, 
W5  =  2.255mm  and  L5  =  4.719mm. 

The  three  director  system  was  also  fabricated  on  0.5mm 
thick  Trans-Tech  D4  material,  and  the  ferrite  material  was 
Trans-Tech  TTI-3000  with  a  radius  of  1.75mm.  Once  again, 
the  externally  applied  bias  field  was  set  at  2575  Oe.  The  match¬ 
ing  networks  connected  to  the  two  50Q  ports  are  identical 
to  those  of  the  two-director  system.  The  matching  network 
connecting  the  circulator  to  the  three-director  antenna  has 
the  following  geometrical  parameters:  W\  =  1.2mm,  L\  = 
1.101mm,  W2  =  0.216mm,  L2  =  1.427mm,  W3  =  1.561mm, 
L3  =  6.473mm,  W4  =  0.216mm,  L4  =  11.47mm,  Wb  = 
2.255mm  and  L5  =  4.719mm. 

The  measured  and  simulated  data  for  the  final  fabricated 
systems  are  shown  in  Figures  15  and  16.  The  measured  return 
loss  for  each  device  is  about  25  dB;  the  measured  isolation 
varies  from  20  dB  to  35  dB  across  the  band  of  interest.  The 
insertion  loss  is  about  0.8  dB  and  is  quite  flat  across  the 
band.  These  measured  results  correlate  quite  well  to  results  of 
simulation.  Moreover,  we  regard  isolations  and  return  losses 
in  excess  of  20  dB  to  be  quite  good  and  in  excess  of  25  dB 
to  be  superb. 


Fig.  15.  Measured  (solid  lines)  and  simulated  (dashed  lines)  scattering 
parameters  of  the  two-director  element  system  with  a  2575  Oe  applied  bias. 


Although  the  plots  shown  in  Figures  15  and  16  illustrate 
the  results  of  the  as-designed  systems,  improved  performance 
may  be  attained  by  varying  the  applied  bias  field.  For  example. 
Figure  17  shows  the  measured  scattering  parameters  of  the 
two-director  antenna  system  with  an  applied  bias  of  2291 
Oe  rather  than  that  of  the  previous  measurement  of  2575  Oe. 
These  tuned  parameters  exhibit  a  very  good  isolation  metric 
of  30  dB  over  a  bandwidth  of  around  6%.  Similarly,  Figure 
18  shows  the  measured  scattering  parameters  of  the  three- 
director  antenna  system  with  an  applied  bias  of  2280  Oe  rather 
than  that  of  the  previous  measurement  of  2575  Oe.  These 
parameters  show  a  very  good  isolation  metric  of  30  dB  over  a 
relatively  narrow  bandwidth  of  2.1%,  and  a  superb  isolation  of 
40  dB  over  a  narrow  bandwidth  of  0.8%.  Hence  by  varying  the 
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Fig.  16.  Measured  (solid  lines)  and  simulated  (dashed  lines)  scattering  Fig.  18.  Measured  scattering  parameters  of  the  three-director  element  system 

parameters  of  the  three-director  element  system  with  a  2575  Oe  applied  bias.  with  a  2280  Oe  applied  bias. 


applied  bias  field,  we  can  dynamically  optimize  these  systems 
for  a  variety  of  applications. 


Fig.  1 7.  Measured  scattering  parameters  of  the  two-director  element  system 
with  a  2291  Oe  applied  bias. 

The  optimized  bias  point  is  shown  on  the  measured  hystere¬ 
sis  curve  of  Figure  1 3,  which  shows  that  the  material  is  in  a 
partially  saturated  state.  This  partial  saturation  is  significant  in 
light  of  the  commonly  held  belief  that  losses  increase  as  the 
material  moves  out  of  saturation.  The  data  shown  in  Figure 
17  and  18  illustrate  that  this  assumption  does  not  always 
hold,  at  least  for  our  application.  At  15  GHz,  for  example, 
the  insertion  loss  changes  from  0.8  dB  to  0.68  dB  for  the 
two  director  system  and  from  0.72  dB  to  0.68  dB  for  the 
three  director  system.  Hence  losses  have  reduced,  rather  than 
increased,  when  the  ferrite  material  is  partially  saturated. 

To  illustrate  the  repeatability  of  these  designs,  two  circuits 
were  fabricated  of  each  system  design.  The  measurements  of 
each  individual  circuit  showed  very  good  correlation,  which 
points  to  the  insensitivity  of  the  frequency  response  from 
device  to  device. 


Conclusion 

Two  Ku-band  circulator/antenna  systems  were  presented 
which  exhibit  superb  performance  in  the  desired  frequency 
spectrum.  Each  system  was  magnetically  tuned  to  optimize 
isolation  performance.  The  two  director  system  realized  6% 
bandwidth  for  a  30  dB  isolation  metric  with  return  loss  better 
than  20  dB.  The  three  director  system  realized  2.1%  bandwidth 
for  a  30  dB  isolation  metric  and  0.8%  bandwidth  for  a  40  dB 
isolation  metric;  in  both  cases  the  return  loss  remained  better 
than  20  dB.  The  measured  isolation  metric  better  than  40  dB 
allows  for  the  design  of  narrow  band  simultaneous  transmit 
and  receive  systems  in  Ku  band. 
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The  Circulator  and  Antenna  as  a  Single  Integrated 

System 

Ryan  S.  Adams,  Member,  IEEE,  Benton  O’Neil,  Member,  IEEE,  Jeffrey  L.  Young,  Fellow,  IEEE 


Abstract — An  integration  methodology  is  presented  for  antenna 
systems  that  incorporate  a  ferrite  circulator  wherein  the  geo¬ 
metrical  and  material  parameters  of  the  circulator  are  adjusted 
to  effect  integration;  no  impedance  equalizers  are  required.  To 
validate  the  proposed  methodology,  a  system  comprised  of  a 
microstrip  circulator  and  Yagi  planar-type  antenna  is  presented. 
The  overall  system  exhibits  excellent  passband  characteristics 
in  both  the  experimental  and  numerical  data.  Additionally, 
the  ferrite  properties  and  system  network  parameters  were 
extensively  measured  to  substantiate  the  claims  made  herein. 

Index  Terms — Ferrites,  Microwave  Circulators,  Microstrip  An¬ 
tennas. 


I.  Introduction 

When  high-performance  RADAR  or  communication  sys¬ 
tems  are  designed  it  is  common  to  break  the  system  into  its 
major  components  and  design  each  to  an  interface  specifica¬ 
tion;  each  component  can  then  be  addressed  separately  from 
the  rest  of  the  system.  The  interface  specification  for  each 
component  is  typically  chosen  to  be  a  purely  real  impedance 
equal  to  some  standardized  characteristic  impedance  (e.g., 
50  12).  This  approach  simplifies  the  overall  design  process, 
but  it  does  not  assure  optimality  or  efficiency  of  the  overall 
system  and  typically  requires  excessive  real  estate  in  the  final 
implementation. 

In  the  present  work  we  show  that  a  system  that  incorporates 
a  circulator  (such  as  a  monostatic  RADAR  system)  can  be  de¬ 
signed  without  the  need  for  impedance  equalization.  This  task 
is  accomplished  by  the  appropriate  selection  of  geometrical 
and  material  parameters  to  establish  the  proper  characteristic 
impedances  and  circulation  impedances  at  each  port  to  natu¬ 
rally  interface  with  the  adjoining  system  components. 

To  verify  the  approach  herein  a  system  is  designed  that 
consists  of  a  circulator  and  a  planar  Yagi  antenna.  We  adopt 
the  following  notation:  port  1  of  the  circulator  is  designed 
to  interface  with  the  antenna  and  ports  2  and  3  are  designed 
to  interface  with  a  standard  50  12  characteristic  impedance  to 
allow  accurate  measurement  of  system  network  parameters. 
The  measured  system  response  reveals  25  dB  isolation  with 
14  dB  return  loss  and  less  than  2.5  dB  insertion  loss  over  a 
bandwidth  of  10.5%  centered  at  9.5  GHz. 

Manuscript  received  XXXX  XX,  20XX;  revised  XXXX  XX,  20XX.  This 
work  was  supported  by  the  Office  of  Naval  Research  under  Award  Number 
N000 14—04- 1-0272. 

Ryan  S.  Adams  is  with  the  Department  of  Electrical  and  Computer 
Engineering,  University  of  North  Carolina  Charlotte,  Charlotte,  North  Carolina 
(email:  radams41@uncc.edu). 

Benton  O’Neil  and  Jeffrey  L.  Young  are  with  the  Department  of  Electrical 
and  Computer  Engineering,  University  of  Idaho,  Moscow,  Idaho. 


Fig.  1.  Geometrical  parameters  associated  with  the  circulator  design. 

II.  Circulator  Modeling 

The  intent  of  this  work  is  to  exploit  the  material  and 
geometrical  parameters  of  a  circulator  to  create  an  integrated 
system  with  desirable  circulation  and  radiation  characteristics, 
without  the  use  of  impedance  equalizers.  A  circulator  is  a 
three  port  device  that  consists  of  a  cylindrical  ferrite  “puck” 
with  transmission  line  ports  emanating  radially  outward.  The 
geometrical  parameters  associated  with  this  device  arc  radius 
a,  coupling  angle  t/\  thickness  h  and  port  angle  <j>,  as  illustrated 
in  Figure  1.  The  material  parameters  are  magnetic  saturation 
An Ms,  effective  internal  magnetic  intensity  H0 ,  ferrite  dielec¬ 
tric  constant  tf  and  dielectric  constant  of  the  surrounding 
material 

In  traditional  circulator  design  the  coupling  angles  and  port 
angles  are  chosen  to  be  identical  for  each  port,  which  assures 
that  the  load  impedance  for  perfect  circulation,  as  presented  in 
[1],  is  identical  at  each  port.  When  this  is  the  case  the  trans- 
impedance  Green’s  function  outlined  in  [2]  may  be  directly 
employed  to  determine  network  parameters  of  the  device.  To 
integrate  the  circulator  with  an  antenna  the  Green’s  function 
of  [2]  must  be  slightly  modified  to  account  for  non-identical 
port  and  coupling  angles.  For  this  case  the  elements  of  the 
impedance  matrix  are  given  by 

Z{j  =  /  (1) 

where  G  (&,</>')  is  the  Green’s  function  of  [2].  All  nine 
impedance  parameters  of  the  circulator  must  be  computed  to 
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fully  characterize  the  device  since  the  coupling  angles  and  port 
angles  are  not  rotationally  symmetric. 

To  avoid  an  abrupt  transition  in  the  microstrip  trace,  each 
port  of  the  circulator  is  chosen  to  equal  the  width  of  the  feed 
trace  of  its  corresponding  system  component.  For  example,  if 
the  width  of  the  feed  structure  of  the  component  at  port  1 
equals  w\,  then  the  width  of  the  circulator  at  port  1  must  also 
equal  w\  and 

^=8in-l(§)-  (2) 
A  similar  expression  is  used  for  ports  2  and  3.  This  equation 
clearly  shows  that  the  radius  of  the  puck  may  be  varied  to 
effect  system  design,  but  the  coupling  angles  of  the  individual 
ports  may  not. 

The  load  impedances  seen  by  each  port  of  the  circulator  may 
be  complex  and  frequency  dependent,  so  the  interconnecting 
microstrip  lengths  may  contribute  to  a  broadband  response. 
For  this  reason  L\%  L2  and  L3  of  Figure  1  are  included  as 
variables  in  the  design  process. 


III.  Design  Approach 


The  system  integration  procedure  begins  with  a  full  charac¬ 
terization  of  the  components  that  will  connect  to  the  circulator. 
This  characterization  can  be  accomplished  with  a  full-wave 
solver,  measurement,  or  existing  model.  To  assure  a  valid 
characterization,  all  components  must  incorporate  a  common 
substrate  dielectric  constant  and  thickness.  Additionally,  each 
component  must  have  a  well  defined  input  microstrip  trace 
geometry  that  will  remain  constant  throughout  the  process. 

To  complete  the  system  integration  procedure,  the  Green’s 
function  is  employed  as  described  previously  to  choose  appro¬ 
priate  values  of  microstrip  interconnecting  trace  lengths  Li, 
L2,  L3,  radius  a,  port  angle  fa,  fa  and  saturation  magne¬ 
tization  47 rMs;  these  parameters  are  shown  in  the  diagram 
of  Figure  1.  The  Green’s  function  is  invoked  in  connection 
with  an  appropriate  optimization  procedure  to  determine  an 
acceptable  system  response. 

For  this  integration  approach  we  have  chosen  to  employ  a 
multi-objective  optimization  procedure  based  on  the  standard 
circulator  metrics  of  return  loss,  insertion  lass  and  isolation. 
With  the  load  reflection  coefficients  given  by  I\,  t  =  1,2,3, 
the  return  loss  at  port  i  is  given  by 


RU  = 


(1  -  |r\|2)  (1  -  iiwi2) 


|i-rjv«r 


(3) 


and  the  power  delivered  to  port  i  from  port  j  is  given  by 


n~ 

zoj  1 

(1  -  ILI2) 

1(1 -IUI2) 

*5" 

Z0{ 

|1  —  TjTinjl 

(4) 


In  each  of  these  equations,  Tiriii  is  the  input  reflection  co¬ 
efficient  at  port  i  of  the  circulator  when  the  remaining  ports 
are  loaded  with  their  respective  reflection  coefficients,  ZQ{  is 
the  port  characteristic  impedance  of  the  ith  port  and  V*  is  the 
voltage  at  port  i  when  each  port  is  loaded  with  its  respective 
reflection  coefficient. 

It  should  be  noted  that  the  Green’s  function  provides  a  good 
design  tool,  as  it  is  computationally  robust,  but  the  results 


Fig.  2.  Simulated  E-plane  antenna  patterns  of  the  antenna  of  [31  with  various 
substrate  thicknesses. 


are  approximate.  Consequently,  exclusive  use  of  the  Green’s 
function  for  system  design  will  create  a  design  that  operates 
near  the  desired  value,  but  may  not  meet  specifications.  Further 
design  modifications  may  be  required  with  a  full-wave  solver 
prior  to  fabrication. 


IV.  Results 

A  system  comprised  of  a  ferrite  circulator  and  planar 
antenna  was  designed  and  fabricated  following  the  theory 
presented  above.  The  antenna  was  chosen  to  interface  to 
port  1  of  the  circulator;  ports  2  and  3  were  designed  to 
interface  to  a  standard  50  il  characteristic  impedance  to  aid 
in  accurate  network  parameter  measurement.  The  substrate 
material  used  in  this  system  design  was  USTC995  by  US 
Technical  Ceramics,  which  has  a  dielectric  constant  of  fd  = 
9.6  and  tan  S  =  0.0007. 

The  antenna  element  was  chosen  from  among  three  options 
based  on  the  planar  Yagi  antenna  of  Kaneda  [3];  modifica¬ 
tions  were  made  in  the  substrate  thickness,  dielectric  constant 
(ed  =  9.6)  and  minor  changes  in  the  miters.  The  substantive 
difference  between  the  three  antennas  was  the  choice  of 
dielectric  thickness;  values  of  0.5  mm,  0.75  mm  and  1 .0  mm 
were  investigated.  The  antenna  with  each  substrate  thickness 
was  simulated  using  Ansoft’s  HFSS  to  determine  pattern  and 
impedance  data.  Figures  2  and  3  show  the  E-plane  and  H- 
plane  antenna  patterns,  respectively,  and  Figure  4  shows  the 
input  impedance  values  of  each  antenna.  The  trace  width  of 
the  antenna  feed  line  in  all  cases  was  fixed  at  1.2mm,  which 
also  became  the  width  of  port  1  of  the  circulator. 

The  optimum  criteria  for  the  choice  of  antenna  was  a  flat 
impedance  response,  with  approximately  zero  imaginary  part 
near  10  GHz,  as  this  response  was  believed  to  be  easily 
matchable  with  a  circulator  at  that  frequency.  In  addition  to 
this  impedance  criteria,  it  was  desired  that  the  antenna  pattern 
exhibit  a  large  main  beam  with  low  sidelobes  in  both  principal 
planes.  The  antenna  with  0.75mm  dielectric  was  chosen  for 
integration,  since  it  exhibited  better  than  20dB  side  lobe  levels 
and  the  best  impedance  response  of  the  available  choices.  This 
choice  of  dielectric  constant  and  thickness  corresponds  to  a  50 
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Fig.  3.  Simulated  //-plane  antenna  patterns  of  the  antenna  of  [3]  with 
various  substrate  thicknesses. 


Fig.  4.  Simulated  input  impedance  of  the  antenna  of  (3]  with  a  substrate 
thickness  of  0.75mm. 


Q  microstrip  line  width  of  0.742mm,  which  defined  the  widths 
of  ports  2  and  3  of  the  circulator. 

The  resulting  circulator  design  is  shown  in  Figure  5,  and 
the  fabricated  system  is  shown  in  Figure  6.  The  geometrical 
parameters  for  the  circulator  were  found  to  be  a  =  2.6mm,  h 
=  0.75mm,  t/>i  =  0.233  rad,  V>2  =  =  0.14  rad,  02  =  2.15 

rad,  03  =  -2.15  rad  and  Lx  =  2.9mm.  In  this  design,  L2  and 
L3  were  arbitrarily  chosen  since  ports  2  and  3  were  loaded 
with  their  respective  characteristic  impedances.  For  the  test 
circuit  shown  in  Figure  6,  L2  =  £3  —  20mm  as  measured 
from  the  edge  of  the  dielectric  to  the  edge  of  the  puck.  The 
ferrite  material  chosen  for  this  system  was  Trans-Tech  TT1- 
1000  with  material  parameters  47 xMa  -  1000  Gauss,  e/  =  1 1 .6 
and  tan  Sf  =  0.00025. 

The  ferrite  hysteresis  and  ferromagnetic  resonance  proper¬ 
ties  were  measured,  with  results  shown  in  Figures  7  and  8, 
respectively.  The  hysteresis  curve  demonstrates  that  although 
the  material  exhibits  a  saturated  magnetization  of  1000  Gauss, 
the  magnetization  at  the  bias  point  is  approximately  900 
Gauss-a  reduction  of  10%  from  full  saturation.  Additionally, 
the  peak-to-peak  linewidth  of  the  material  was  measured  at 
10.0  GHz  to  be  103  Oe  which  corresponds  to  half  power 
linewidth  of  179  Oe.  As  will  be  seen  in  the  ensuing  simulation 


Port  1 


Fig.  5.  Puck  design  used  with  the  0.75mm  thickness  antenna  connected  to 
port  1. 


Fig.  6.  Overall  system  comprised  of  the  feirite  puck  of  Figure  5  and  the 
antenna  of  [3]  with  a  0.75mm  thick  dielectric. 


data  this  linewidth  value  leads  to  moderately  low  insertion 
losses. 

Comparison  of  the  simulated  and  measured  response  is 
shown  in  Figures  9  and  10.  These  figures  clearly  illustrate  a 
circulating  response,  with  a  return  loss  above  14  dB,  insertion 
loss  below  2.5  dB  and  superb  isolation  better  than  25  dB  over 
the  band  from  9-10  GHz.  Certain  deviations  from  simulation 
are  apparent  particularly  in  the  insertion  loss  measurement 
which  was  simulated  to  be  in  the  vicinity  of  0.5  dB.  First, 
the  higher  than  expected  insertion  loss  is  linked  to  the  coax 
to  microstrip  transitions  and  fabrication  deficiencies.  Simu¬ 
lation  of  the  circuit  with  transitions  versus  without  showed 
approximately  ldB  insertion  loss  directly  associated  with  the 
use  of  these  transitions.  Second,  the  porosity  of  the  alumina 
substrate  sample  required  that  a  very  thin  metal  layer  be  used 
to  form  the  transmission  lines.  Third,  small  gaps  existed  at  the 
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Fig.  7.  Measured  hysteresis  curve  of  TT1 -1000. 


ferrite/dielectric  boundary  which,  coupled  with  the  thin  metal 
layer,  created  small  gaps  in  the  feed  traces  at  the  boundary  of 
the  circulator.  Fourth,  the  quasi -saturated  state  of  the  ferrite 
per  Figure  7  is  known  to  affect  loss.  Our  model  assumes  a 
fully  saturated  ferrite. 

V.  Conclusion 

A  novel  design  procedure  has  been  devised  for  antenna 
systems  that  incorporate  a  ferrite  circulator.  The  procedure 


Fig.  10.  Simulation  and  measurement  insertion  loss  and  isolation  of  the 
system  of  Figure  6. 


utilizes  the  inherent  material  and  geometrical  parameters  of 
the  circulator,  combined  with  a  Oreen  s  function  model,  search 
methodology  and  appropriate  search  metrics.  This  approach 
significantly  reduced  the  amount  of  real  estate  required  for 
the  integrated  device,  and  reduced  the  copper  losses  associated 
with  the  use  of  open-circuit  stub  based  impedance  equalizers. 

We  have  shown  in  the  present  work  that  impedance  equal¬ 
ization  may  be  accomplished  by  direct  modification  of  the 
circulator  material  and  geometrical  parameters.  One  question 
that  arises  in  connection  with  this  work  is  whether  other 
functionality  can  be  accomplished  by  the  circulator  directly. 
For  example,  the  circulator  is  shaped  and  fed  similar  to  a 
circular  microstrip  patch  antenna;  can  it  be  made  to  radiate 
in  addition  to  circulate  and  provide  impedance  matching  to 
a  system?  This  possibility  would  enable  significant  reductions 
in  real  estate  of  phased  array  antennas  for  monostatic  RADAR 
applications. 
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A  Generalized  Time-Stepping  Algorithm  for 
Electromagnetic  Wave  Propagation  in  Ferrite 
Media 

Ryan  S.  Adams  and  Jeffrey  L.  Young 


Abstract—  A  time  integration  scheme  is  considered  herein  for 
simulations  of  electromagnetic  waves  within  a  fully  saturated  fer¬ 
rite  medium.  By  considering  Maxwell’s  equations  in  conjunction 
with  the  torque  equation  of  motion  for  electron  precession,  we 
show  that  a  stable  second-order  scheme  can  be  achieved  using 
both  central  differences  and  averages.  The  scheme  is  validated 
using  a  stripline  phase  shifter  as  a  test  case.  Not  only  does  the 
simulation  agree  quite  well  with  results  produced  by  a  popular 
frequency-domain,  finite  element  solver,  the  data  also  exhibits  a 
very  desirable,  low-loss,  phase  shift  effect. 

Index  Terms—  Ferrites,  Ferrite  Phase  Shifters,  Magnetic  Ma¬ 
terials,  FDTD. 


I.  Introduction 

During  the  past  decade  there  has  been  considerable  research 
activity  in  the  modeling  of  electromagnetic  waves  within  a 
fully  saturated  ferrite  using  time-domain,  Maxwellian  field 
solvers.  Two  popular  approaches  have  been  considered  to 
date  -  recursive  convolution  (RC)  and  direct  integration  (DI). 
The  RC  methods  use  the  time-domain  representation  of  the 
Polder  terms  of  the  susceptibility  matrix  in  conjunction  with 
the  constitutive  convolution  integral  as  the  governing  equation 
for  the  ferrite.  Depending  on  whether  the  convolution  integral 
is  discretized  using  piecewise  constant  or  piecewise  linear 
interpolation,  Schuster  and  Luebbers  [1]  and  others  were  able 
to  achieve  stable  schemes  that  yielded  first  or  second  order 
solutions.  Other  researchers,  such  as  [2]  -  [4],  adopted  vari¬ 
ous  direct  integration  methods,  whereby  the  dynamic  torque 
equation  of  electron  precession  is  used  as  the  starting  point.  In 
these  cases,  one-sided  differences,  central  differences  and/or 
central  averages  are  invoked  to  discretize  the  torque  equation. 

The  scheme  considered  herein  continues  the  development 
of  the  direct  integration  methodologies  by  following  a  similar 
line  of  thought  advocated  by  [4].  However,  instead  of  rewriting 
the  torque  equation  in  terms  of  the  magnetic  flux  density  B 
and  the  magnetic  intensity  H,  as  was  done  by  [4],  we  opt 
to  write  the  torque  equation  in  terms  of  the  magnetic  flux 
density  B  and  the  magnetization  M.  By  doing  so,  we  leave 
the  torque  equation  in  its  natural  state  variable  form,  which 
allows  us  then  to  discretize  it  using  central  differences  and 
averages.  The  resulting  discretized  equation  is  couched  in  a 
compact  form  using  dyadic  and  vector  operations  and  is  valid 
for  any  direction  of  ferrite  bias  field. 

It  has  been  our  goal  to  follow  the  philosophy  of  developing 
FDTD  schemes  for  dispersive  media  by  adding  equivalent 
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magnetization  terms  to  Faraday’s  law  (as  in  the  case  of  femtes) 
or  polarization  terms  to  Ampere’s  law  (as  in  the  case  of  Debye 
[5],  Lorentz  [6]  or  plasma  media  [7]).  By  doing  so,  we  never 
modify  the  core  Yee  algorithm,  as  expressed  in  terms  of  E  and 
B  in  Faraday’s  law  and  H  and  D  in  Ampere’s  law,  regardless 
of  the  media  type  considered.  The  scheme  described  herein 
achieves  this  goal. 

To  validate  the  proposed  integration  methodology,  a  three- 
dimensional,  stripline  phase  shifter  is  considered  in  which  the 
field  bias  is  in  the  direction  of  propagation.  This  device  is 
particularly  attractive  for  our  purposes,  since  various  biasing 
fields  need  to  be  considered  to  see  the  phase  shift  effect.  Using 
a  solution  produced  by  a  frequency-domain,  finite-element 
solver  as  a  standard,  we  show  near  perfect  agreement  between 
that  standard  and  the  data  produced  by  the  proposed  scheme. 
Moreover,  from  a  device  performance  point  of  view,  a  full 
180°  phase  shift  can  be  observed  by  adjusting  the  effective 
internal  field  from  0  Oe  to  2800  Oe;  the  return  loss  remains 
below  20  dB  and  the  transmission  loss  remains  below  0.5  dB 
over  this  same  range.  To  validate  that  the  proposed  scheme 
can  accomodate  an  arbitrary  direction  of  magnetic  bias,  this 
geometry  was  again  simulated  using  an  obliquely  directed 
bias  direction;  excellent  correlation  is  again  shown.  As  a  final 
analysis  of  this  device,  a  simulation  was  completed  using 
1,000,000  timesteps;  no  late-time  stability  problems  were 
noted. 


II.  Formulation 


Within  the  ferrite  material,  the  magnetization  M  and  mag¬ 
netic  intensity  H  are  related  to  each  other  by  the  Gilbert 
equation  of  motion,  including  the  phenomenological  loss  term 
[8].  This  equation  can  be  expressed  as 


dM 

dt 


a  _  _  dM 

=  x  H  -  x 


(1) 


where  /x0  is  the  permeability  of  free  space,  7  =  1.759  x  101 1 
C/kg  is  the  gyromagnetic  ratio  and  a  is  a  damping  factor. 
The  vectors  M  and  H  consist  of  both  dc  and  ac  components. 
The  dc  component  of  H  is  the  total  effective  internal  dc 
magnetic  intensity;  the  dc  component  of  M  is  the  internal 
dc  magnetization. 

The  ac  magnetic  flux  density  is  related  to  the  magnetic 
intensity  and  magnetization  as  [9] 

Har  =  — -B„c  —  M„r,  (2) 

Ho 

so  the  total  magnetic  intensity  within  the  ferrite  is  given  by 

H  =  +  —  Bac  —  M0(V  (3) 

Mo 

In  this  work  we  define  a  vector  Larmor  precession  frequency 
cJ0  and  vector  magnetization  frequency  u;m  such  that 


uQ  =  /*o7H  dc  =  v0xa*  +  Uoy&y  +  Uot*z  (4) 


and 


UJm  =  Mo7Mdc  =  CJmxaT  + 


(5) 
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Eqn.  (1)  can  then  be  expressed  as 

0Mac 


and 


dt 


—  —  fJiolNldc  X  H dc  -  Mac  X  (^rn  +  <*>o) 


+  — Boc 
Mo 

a  <9M 

+  = 


X  LJm  —  7 Mac  X  Ba 


Wfl 


dt 


x  u,n  - 


a 


M„ 


0Ma 


(6) 


iMr*ac~  ot 

If  we  assume  that  |MaC|  <  |MdC|,  |Hac|  -C  |HdC|,  M^c  is 
parallel  with  Hac  and  |M<fc|  —  M9,  then  Eqn.  (6)  becomes 

,  _  a  ^MflC  __ 

t  h^acX(cJtn  "1“  ^o)“H  , 

"  (7) 


^Mac  ^  d  xrr 
—  - JL>ac  X  Ll) 

Mo 


dt 


In  so  far  as  the  equation  of  motion  given  by  Eqn.  (1)  is 
valid,  Eqn.  (7)  represents  the  relationship  between  M  and  B 
within  a  fully  saturated  ferrite  under  small-signal  conditions. 
This  equation,  as  cast  in  terms  of  dc  biasing  fields  and  ac 
signal  fields  will  be  the  basis  for  the  subsequent  discretization 
procedure.  For  the  purposes  herein,  we  will  regard  the  Gilbert 
form  of  the  Landau-Lifschitz  model  to  be  phenomenologically 
true  in  the  form  given  by  Eqn.  (7)  [10]. 

Along  with  Eqn.  (7),  the  equations  governing  the  electro¬ 
magnetic  problem  are: 

0Eqc  (8) 


dt 

dBoc 

dt 


=-V 

6 


=  -  V  X  Eac 


(9) 


and  Eqn  (2),  where  Eac  is  the  ac  electric  intensity  and  Bac  is 
the  ac  magnetic  flux  density.  These  equations,  along  with  Eqn. 
(7)  can  be  temporally  discretized  using  both  central  differences 
and  averages  as  follows: 

E"+1/2  =  Kcl/2  +  -  V  x  H”c, 


B”c+1=B”c-5(VxE2c+1/2, 

+  *t_  /BS+1  +  B 


(10) 

(ID 
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N 

£ 
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1 

M2  = 


1  M2,  —  M2, 

-M2.  1  m2. 

Mo  -M2.  1 


(16) 


The  elements  of  Mi  and  M2  are  given  by 

,n  =  x,y,z  (17) 


«  ^  fit  /  .  \  auJmn 

MXn  =  T  (Wm'‘+WoJ“pJ 


and 


M2n  =  --(um„+“oJ  ,_n 


,n  =  x,y,  z.  (18) 


With  these  definitions,  Eqn.  (12)  can  be  expressed  as 

£-%'  sm ■  (b;*1  +  , 

^  (19) 

where,  for  example,  M2  •  M£c  represents  a  dyadic-vector  dot 
product  operation. 

To  complete  our  treatment,  we  choose  to  employ  the  Yee 
cell  of  Figure  1  in  conjunction  with  spatial  central  differences 
and  averages.  With  this  cell  arrangement,  Eqns.  (10)  -  (1 1)  can 
be  spatially  discretized  as  in  the  traditional  Yee  scheme  [11]. 
Since  the  components  of  the  magnetization  vector  have  been 
placed  coincident  with  the  magnetic  flux  density  and  magnetic 
intensity,  Eqn.  (13)  requires  no  spatial  averaging.  With  the  grid 
chosen  according  to  Figure  1,  Eqn.  (19)  shows  that  M”1 1 
depends  on  the  previous  values  of  M",  M",  By4  and 
B?+1/2,  each  of  which  is  not  located  at  the  same  point  in 
the  grid  as  Af”+1.  Hence,  each  of  these  components  require 
four-point  spatial  averaging  to  compute  the  correct  value  at  the 
location  of  Mx  within  the  grid.  Similar  statements  can  be  made 
about  four-point  spatial  averaging  for  the  update  equations  for 
M”+1  and  M"+1. 


(12) 

(13) 


H"r+1  =  -BS+1  -  M"r+1. 

Mo 

It  should  be  noted  that  Moc,  Hoc  and  Bar  are  all  temporally 
synchronized  in  the  aforementioned  scheme.  Equations  (10) 
-  (13)  are  valid  for  any  numerical  scheme  that  employs  both 
temporal  central  differences  and  averages. 

To  obtain  an  implementation  algorithm  associated  with  Eqn. 
(12),  we  must  first  solve  for  M”c+1.  To  do  this,  we  define  the 
following  dyads  couched  in  matrix  form: 


(14) 


(15) 


III.  RESULTS 

Consider  the  stripline  phase  shifter  topology  shown  in 
Figure  2.  This  problem  is  an  appealing  validation  geometry 
for  the  algorithm  presented  herein,  as  it  provides  a  mixture 
of  different  materials  (substrate,  PEC  and  ferrite),  allows  a 
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comparison  of  several  different  values  of  effective  internal  field 
and  requires  an  applied  bias  in  a  direction  different  from  bias 
normal  to  the  plane  of  the  circuit. 

The  geometrical  parameters  chosen  for  this  validation  are 
Wf  =  2.9mm,  Lf  =  24mm,  W  =  1.1mm  and  h  =  2mm;  the 
width  of  the  domain  is  8.6mm.  The  ferrite  material  chosen 
is  Trans-Tech  TTVG-1000  which  has  a  dielectric  constant 
of  6/  =  14,  magnetic  saturation  AttM8  —  1000  G  and 
linewidth  A H  =  10  Oe,  as  measured  at  fmeas  =  9.4GHz.  The 
surrounding  material  is  Trans-Tech  D-16  Mg-Ti  which  has  a 
dielectric  constant  of  ej  —  16.  From  these  parameters,  ujm  = 
li0~fAds  —  1.76  x  1010  rad/s  and  a  =  n0'yAH/2umeaa  = 
0.237  x  10-4.  The  width  of  the  trace  is  set  to  achieve  a  250 
characteristic  impedance. 


wavelength  is  0.158mm.  For  convenience  in  accomodating  the 
correct  dimensions  described  previously,  the  following  cell 
dimensions  were  chosen:  6X  =  0.15mm,  6y  —  0.1mm  and 
5Z  =  0.1mm.  The  overall  domain  size  was  chosen  to  be  100 
cells  in  the  y-< iirection,  180  cells  in  the  s -direction  and  20 
cells  in  the  z-direction.  The  timestep  was  chosen  utilizing  a 
CFL  stability  criteria  of  0.5  so  that  St  «  0.1  ps. 

To  validate  the  algorithm  presented  herein,  a  comparison 
of  the  real  and  imaginary  transmitted  scattering  parameters 
of  the  phase  shifter  was  made  with  Ansoft’s  High  Frequency 
Structure  Simulator  (HFSS),  as  shown  in  Figure  3.  In  this  plot, 
we  see  very  good  correlation  between  the  two  data  sets,  thus 
indicating  that  for  the  fully  saturated  and  low  power  case,  the 
algorithm  is  consistent  with  frequency-domain  formulations. 


Fig.  2.  Stripline  phase  shifter  validation  topology.  The  darker  shaded  region 
represents  the  ferrite. 


The  stripline’s  TEM  waves  are  excited  and  detected  per 
the  method  described  in  [12].  This  approach  allows  for  a 
much  smaller  domain  size  than  would  otherwise  be  required 
as  the  source  and  terminal  planes  of  the  simulation  can  be 
placed  arbitrarily  near  the  ferrite/dielectric  interface  without 
loss  of  accuracy.  The  pulse  chosen  to  characterize  the  circuit 
is  gaussian  of  width  a,  where  a  =  6.065  x  1010  l/s;  the  pulse 
is  of  the  form  e_<2//a2. 

To  determine  appropriate  cell  size,  we  must  first  consider 
the  wavelength  within  the  ferrite  material.  Per  [10],  A  =  2n/k 
and  k  =  where  /ie  =  -  «2)//*/,  and 


=  H  o 


MpLdfn  \ 

u2-u2)' 


K  =  Ho 


u2q-u2' 


(20) 

(21) 


At  an  operating  frequency  of  13  GHz,  A  =  6.3mm  for  an 
effective  internal  field  of  zero,  and  A  =  8.1mm  for  an  effective 
internal  field  of  3000  Oe.  To  achieve  superb  accuracy,  the  cell 
size  is  chosen  to  be  less  than  A/40  which  for  the  shortest 


Fig.  3.  A  comparison  of  S21  with  HFSS  for  effective  internal  field  equal  to 
zero.  Solid  lines  are  FDTD  and  dashed  lines  are  HFSS. 

The  simulated  return  loss  of  this  device  at  13  GHz  is  shown 
in  Figure  4.  This  plot  indicates  that  the  device  is  usable  in  a 
250  system  with  better  than  20dB  specification  for  an  effective 
internal  field  less  than  2900  Oe.  The  simulated  phase  shift  at 
13  GHz  is  shown  in  Figure  5  and  shows  that  for  a  reasonable 
effective  internal  field  strength,  this  device  can  realize  greater 
than  180°  phase  shift. 


Fig.  4.  Return  loss  of  the  device  as  a  function  of  applied  field  when  /  =  13 

GHz. 
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Fig.  5.  Phase  shift  of  the  device  as  a  function  of  applied  field  when  /  =  13 
GHz. 

To  achieve  variation  in  the  effective  internal  field  as  shown 
in  Figs  4  and  5,  variation  of  the  externally  applied  bias  is 
necessary.  To  determine  what  the  appropriate  external  field 
should  be,  we  assume  a  homogeneous  ferrite,  in  which  case, 
Hint  =  Happlied  -  (4nMaNyy),  where  Nyy  is  the  demag¬ 
netization  factor.  Due  to  the  geometrical  shape  of  the  ferrite 
material,  the  demagnetizing  factor,  as  given  by  Joseph  and 
Schlomann  [13],  is  approximately  zero,  hence  the  effective 
internal  field  is  equal  to  the  externally  applied  bias. 

One  of  the  significant  contributions  of  the  proposed  scheme 
is  its  ability  to  model  ferrites  whose  bias  direction  points 
in  a  direction  different  from  the  cartesian  coordinate  axes. 
To  validate  this  capability,  the  simulation  parameters  above 
were  used  with  a  ferrite  having  an  oblique  bias  direction  and 
an  internal  field  of  zero.  The  bias  was  chosen  to  have  the 
following  direction  cosines:  cos  a  =  0.54,  cos/?  =  0.31  and 
cos  7  =  0.78,  where  a  is  measured  from  the  x-axis,  (5  from 
the  y-axis  and  7  from  the  2-axis.  A  comparison  of  the  real 
and  imaginary  transmitted  scattering  parameters  of  the  device 
as  determined  from  the  proposed  scheme  and  HFSS  is  shown 
in  Figure  6.  As  before,  this  plot  shows  very  good  correlation 
between  the  two  data  sets.  As  an  additional  validation  of 
the  proposed  scheme,  this  simulation  was  accomplished  with 
1,000,000  timesteps  to  examine  late-time  stability;  no  late¬ 
time  stability  problems  were  noted. 

IV.  Conclusion 

A  highly  efficient  time-domain  scheme  has  been  presented 
for  the  simulation  of  electromagnetic  waves  within  ferrite 
media.  In  this  investigation,  we  limited  our  treatment  to 
electromagnetic  modes  of  low  wavenumbers.  However,  for 
very  large  wavenumbers  (i.e.  k  >  10  r/cm),  it  is  well  known 
that  electric  intensity  becomes  small  enough  to  be  neglected. 
The  magnetic  intensity,  however,  becomes  irrotational  and  the 
mode  that  propagates  is  regarded  as  magnetostatic.  For  even 
larger  wavenumbers  (i.e.,  k  >  105  r/cm),  a  spin  wave  can  be 
supported  in  the  ferrite.  Due  to  these  large  wavenumbers,  it  is 
questionable  whether  the  explicit  scheme  presented  herein  is 
satisfactory,  due  to  the  limitations  of  the  CFL  number.  (This 
problem  was  also  noted  in  [7]  with  respect  to  the  simulation 


Fig.  6.  A  comparison  of  52 1  with  HFSS  for  an  oblique  bias  direction  and 
effective  internal  field  equal  to  zero.  Solid  lines  are  FDTD  and  dashed  lines 
are  HFSS. 

of  waves  in  warm  plasmas.  In  that  case,  the  large  disparities  in 
wavenumbers  between  the  optical  and  acoustic  modes  created 
algorithmic  difficulties.)  Future  work  in  this  area  is  clearly 
needed. 
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